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Fig. 1. Radial profiles across an initially 2-km-high and 150-km-radius 
plateau. The plateau and upper 2.5 km of the crust are of diabase composition 
and the lower 2.5 km of crust is of gabbro and overlies a peridotite mantle. The 
lithosphere is 33 km thick. The thermal gradient through the upper crust is 
28 K km -1 . Topographic profiles are shown for the original plateau, the 
instantaneous elastic response to loading (isostatic balance), and for 0.5. 10, 
50, and 500 Ma. Vertical exaggeration is x21. 


from those of other workers in that (1) our thermal models are 
distinctive and we believe more realistic and (2) our models are 
dynamically isostatically balanced with both hydrostatic restoring 
forces and dynamic support by tectonic stress and basal tractions 
being modeled. The implications of our models * stress evolution for 
surface deformation can be tested against Magellan imagery. 
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FLEXURAL ANALYSIS OF UPLIFTED RIFT FLANKS ON 
VENUS, Susan A. Evans, Mark Simons, and Sean C. Solomon, 
Department of Earth, Atmospheric, and Planetary Sciences, Massa- 
chusetts Institute of Technology, Cambridge MA 02139, USA. 

Introduction: Knowledge of the thermal structure of a planet 

is vital to a thorough understanding of its general scheme of 
tectonics. Since no direct measurements of heat flow or thermal 
gradient are available for Venus, most estimates have been derived 
from theoretical considerations or by analogy with the Earth [1] . 
The flexural response of the lithosphere to applied loads is sensitive 
to regional thermal structure. Under the assumption that the yield 
strength as a function of depth can be specified, the temperature 


gradient can be inferred from the effective elastic plate thickness 

[2] . Previous estimates of the effective elastic plate thickness on 
Venus range from 1 1—1 8 km for the foredeep north of Uorsar Rupes 

[3] to 30-60 km for the annular troughs around several coronac 
[4,5] . Thermal gradients inferred for these regions are 14-23 K km-* 
and 4-9 K km -1 respectively [3,4]. In this study, we apply the same 
techniques to investigate the uplifted flanks of an extensional rift. 
Hypotheses for the origin of uplifted rift flanks on Earth include 
lateral transport of heat from the center of the rift, vertical transport 
of heat by small-scale convection, differential thinning of the 
lithosphere, dynamical uplift, and isostatic response to mechanical 
unloading of the lithosphere [6]. The last hypothesis is considered 
the dominant contributor to terrestrial rift flanks lacking evidence 
for volcanic activity, particularly for rift structures that are no 
longer active [6]. In this study, we model the uplifted flanks of a 
venusian rift as the flexural response to a vertical end load. 

Tectonic Environment: We examine a linear rift system 
centered at 33°S, 92°E, in an area to the east of Aino Planitia. The 
feature appears as a linear ridge in Pioneer Venus altimetry and thus 
has been named Juno Dorsum. However, the increased resolution of 
Magellan images and topography has established that this feature is 
actually a linear rift with pronounced flanking highs. The rift is 
100 km wide and 450 km long; it has a central depression 1—2 km 
deep and flanks elevated by as much as 1 km (Figs. 1 and 2). The rift 
connects a 1 .5-km-high volcano on its western edge to two coronae, 
Tai Shan and Gefjun [7], to the east. Despite the presence of the 
volcano and coronae, both the center of the rift and the adjacent 
flanks appear lobe free of volcanic flows. Juno Dorsum is at the end 
of a nearly continuous chain of coronae, rifts, and linear fractures 
that extends eastward and appears to terminate at the northwestern 
edge of Artemis Corona, The rift appears in Magellan images as 
numerous east-northeast-trending lineaments, which we interpret to 
be normal faults. To the north and south are smooth, radar-dark 
plains that stand very close to mean planetary elevation. 

We use Magellan altimetric profiles from orbits 965-968 and 
970-972 (Figs. 1 and 2). While topographic profiles perpendicular 
to the strike of the rift are preferable, the orbit tracks cross the rift 



Fig. 1. Topographic contour map of the Juno Donum region. Thin solid 
lines denote positive elevation contours, thick solid lines zero elevation, and 
dashed lines negative elevations. The contour interval is 0.3 km, and the datum 
is mean planetary radius, 6051.9 km. The rift is located in the center right of 
the figure, with an unnamed volcano to its wesL The north- south- trending 
lines indicate the tracks of orbits 965-968 and 971-972 (numbers increasing, 
eastward). The orbit tracks make an angle of approximately 75° with the strike 
of the rift. 
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Fie 2. Example orbit profiles over Juno Dorsum; see Fig-1 for location. 
The rift flank, ate elevated 800 to 1000 m above the surrounding plants and 
interference from other regional topographic features mutmtal. Orb,. 967 
has been offset vertically by 2500 m for clarity. 



at an angle of approximately 75°. and we can easUy make Mte 
geometrical correction for this small difference. Profiles from 
single orbits avoid errors introduced by interorbit offsets. As is 
evident in Fig. 2. the topography of the surrounding plains region is 
rather smooth, and the uplift of the flanks is quite distinct and easy 

to discern from other regional topography. 

Model: We model the uplift of the rift flanks as the flexural 
response of a thin elastic plate overlying a fluid substratum [ 8] . e 
assume that the flexure is of the form 


w = exp 




a) 


where w is the vertical deflection of the plate from the regional 
datum, x is the distance along the profile, and a is the flexural 
parameter given by 
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where E is Young’s Modulus (60 OP.). T e is the elastic plate 
thickness. p m is the mande density (3300 kg uf). g » «£«£«£ 
tional acceleration (8.87 m s 2 ), and v is Poisson s ratio (025). The 
objective is to find the three parameters A, B, and a such that the 
root-mean-square (rms) misfit between the observed and ealeu ated 
topographic profiles is a minimum. Since w depends nonl.nearly on 
a we first fix the value of a, then calculate the best-fitting values 
of A and B. The flexural parameter is then varied over a range 
corresponding to an elastic plate thickness of 1 to 50 km, and the 
procedure is repeated. From the parameters providing the best 
overall fit we calculate the plate curvature at the first zero-crossing^ 
This curvature and the apparent elastic plate thickness are then used 
to estimate the thermal gradient from the curves in [3J. 

We model only the rift flanks, so we remove the centra] nit 
depression from the profiles and model each rift flank separately. 
The location of highest topography in the half profile defines the 
origin of the horizontal axis. We fit topography up to 350 km from 
the rift, to avoid modeling the small rise 400 km to the south of the 
rift (Figs 1 and 2). We also remove the mean elevadon and the best- 
fit linear trend from our profiles before fitting to prevent bias from 

long-wavelength regional topography. 

Results: Of the 14 profiles obtained in orbits 965-968 and 
971-972 we analyze 10 that have flanks rising 600-1200 m above 
the regional mean elevation. Four of these profiles are well fit by 

elastic plate thicknesses of 8-1 6 km. For each of these profiles, there 



a. 3. Top: Topographic prefiles (solid) and best-fit flexural profiles 
ashed) for three half profiles; lop to bottom are otbit 967, southern flank 
bit 968. southern flank orbit 968. nonhem flank Bottom: the nms error. s 
function of elastic plate thickness for the tune profiles. Note that the top two 
rofiles are best fit only over a small range of elastic plate 
re bottom profile cut be fit by a broader range of plate du<^«s« A 
refiles are consistent with elastic plate thicknesses m the ^ 

laritv.the upper two profile, at top have been vertically offset by 1000 m 


a narrow range of plate thicknesses for which the ims error is 
nail, so the results are well constrained (Fig. 3). Three of these 
rofiles produce best-fitting plate thicknesses of 8-9 km. suggest- 
tg that the true effective elastic plate thickness may be near the low 
nd of this range. The remaining profiles generally contain more 
mall-scale topographic variadon of ±100 m in amplitude and are 
qually well fit by a range of plate thicknesses. The fits for all 


8-20 km. 

These effective plate thicknesses and the calculated curvatures 
imply thermal gradients of 15-30 K km ' if the base of the 
mechanical lithosphere is limited by the ducdle strength of ohvme 

[ 31 . Because the values ofelastic lithosphere thickness denv ere 

are comparable to or less than the estimated crustal thickness 
beneath plains regions on Venus [9.10], the base of the mechanical 
lithosphere may lie within the crust and be limited by the ow aw 


would be less than the figures given above. 

The effective elasuc plate thicknesses we derive for the Junonft 
structure (8-20 km) are comparable to those obtained for the North 
Polar Plains at Uorsar Rapes (1 1-18 km) and both regions involve 
flexure of the lithosphere beneath plains at elevations near the 
planetary mean. A general relationship between regional elevation 
and elastic plate thickness might be expected if regional elevauon 
were controlled by lithospheric thermal structure [11], as m the 
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Barth’s oceanic regions. Further analysis of flexure is warranted to 
determine if such a relationship holds for Venus. 
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UNDERSTANDING THE VARIATION IN THE MILLIME- 
TER-WAVE EMISSION OF VENUS. Antoine K. Fahd and 
Paul G. Steffes, School of Electrical Engineering, Georgia Institute 
of Technology, Atlanta GA 30332, USA. 

Recent observations of the millimeter- wave emission from 
Venus at 1 1 2 GHz (2.6 mm) have shown significant variations in the 
continuum flux emission [1 ] that may be attributed to the variability 
in the abundances of absorbing constituents in the Venus atmo- 
sphere. Such constituents include gaseous H 2 S0 4 , S0 2 , and liquid 
sulfuric acid (cloud condensates). Recently, Fahd and Steffes [2,3] 
have shown that the effects of liquid H 2 S0 4 and gaseous S0 2 cannot 
completely account for this measured variability in the millimeter- 
wave emission of Venus. Thus, it is necessary to study the effect of 
gaseous H 2 S0 4 on the millimeter-wave emission of Venus. This 
requires knowledge of the millimeter-wavelength (MMW) opacity 
of gaseous H 2 S0 4 , which unfortunately has never been determined 
for Venus-like conditions. 


We have measured the opacity of gaseous H 2 S0 4 in a C0 2 
atmosphere at 550, 570, and 590 K, at 1 and 2 atm total pressure, and 
at a frequency of 94.1 GHz. Our results, in addition to previous 
centimeter- wavelength results [4], are used to verify a modeling 
formalism for calculating the expected opacity of this gaseous 
mixture at other frequencies. This formalism is incorporated into a 
radiative transfer model to study the effect of gaseous H 2 S0 4 on 
the MMW emission of Venus. 

Experimental Configuration: The experimental setup used 
to measure the MMW opacity of gaseous H 2 S0 4 atmosphere 
consists of a frec-space transmission system as shown in Fig. 1. In 
this system, a glass cell contains the H 2 S 04 /C 0 2 gaseous mixture 
that is introduced prior to the measurement process. The glass cell 
is located inside a temperature-controlled chamber. A transmitting 
antenna is used to launch energy into the glass chamber. A receiving 
antenna is placed at the output of the glass cell in order to collect the 
outgoing signal. Using a precision variable attenuator, the resulting 
opacity of the gaseous mixture is measured. 

Measurement Results: The measured absorptivity (dB/km) of 
^ 2^4 at 94» 1 GHz is shown in Fig. 2 where it is plotted as a function 
of temperature for 2 and 1 atm. The reported absorptive ties in Fig. 2 
are normalized to their respective mixing ratios. The measurements 
were performed at 550, 570, and 590 K in order to allow enough 
H 2 S0 4 vapor pressure in the glass cell. 

Although the measurements were performed at 94.1 GHz, care 
must be taken when projecting the absorption of H 2 S0 4 at frequen- 
cies far from 94. 1 GHz. As a result, we have developed an absorption 
model based on a Van Vleck-Weisskopf (WW) formalism. In this 
formalism, we added the contributions from 2359 resonant lines of 
H 2 S0 4 computed by Pickett et al. (private communication, 1991) 
that cover frequencies between 1 .5 and 450 GHz. 

In order to fully implement the VVW formalism, an appropriate 
broadening parameter must be determined. To solve this problem. 
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Fig. 1 . Block diagram of the atmospheric simulator as configured for measurements of the millimctcr-wavc absorption at 94.1 GHz. 












